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MEASUREMENT  OF  LUNG  FUNCTION  USING  MAGNETOMETERS. 
II.  Data  Acquisition  and  Analysis. 

Donald  L.  Vawter,  Ph.D. 


The  Hyperbaric  Medicine  and  Physiology  researchers  who  measure  lung 
function  using  magnetometers  currently  use  the  program  DATLOG  to  sample  and 
store  the  magnetometer  and  spirometer  output.  To  better  visualize  the 
chest  and  abdomen  during  the  motion  of  respiration,  we  have  developed  a 
software  package  for  plotting  the  sample  data  and  for  fitting  a  variety 
of  models  to  the  data.  There  are  five  basic  programs  discussed  briefly  in 
this  report.  The  programs  are  interactive  in  nature  and  largely  self- 
explanatory.  The  use  of  the  programs  on  one  particular  data  set  is 
presented  in  Appendix  1,  which  should  serve  as  a  user's  guide. 

PROGRAM  ALLDAT 


This  program  prints  the  data  stored  on  the  specified  input  file. 
The  input  file  should  be  in  the  form  generated  by  DATLOG.  The  user 
has  the  option  of  printing  out  data  at  only  a  few  points.  This  program 
is  generally  only  useful  for  inspecting  data  for  spurious  points.  The 
FORTRAN  listing  is  contained  in  Appendix  2. 

PROGRAM  PLT 


Program  PLT  has  three  functions: 

1)  Prepare  a  data  file  for  plotting  the  data  automatically, 
using  program  [1,100]  PLOT.  The  spirometer  output  and  the 
output  from  each  magnetometer  will  be  plotted.  The  user  must 
input  the  name  for  the  plotfile  to  be  generated. 

2)  "DC  Offset"  the  data.  It  may  be  useful  in  some  situations 
to  force  the  data  set  to  be  zero  at  either  the  start  or  the 

end  of  the  data  set.  For  example,  because  the  change  rather  than 
absolute  volume  is  Important,  one  may  wish  to  subtract  the 
spirometer  reading  at  the  first  data  point  from  all  values 
in  the  data  set.  This  is  especially  useful  if  the  first  data 
point  corresponds  to  end-expiratory  values.  The  use  of 
software  to  accomplish  this  function  has  an  advantage  over 
the  previous  method  of  offsetting  the  data  electronically  by 
use  of  a  low-pass  filter  system.  The  advantage  is  that  the 
"raw"  data  is  not  altered.  Hence,  if  one  discovers  later 
that  the  initial  point  does  not  correspond  to  end  expiration, 
one  still  has  the  original  data  set.  The  offset  feature  is 
also  useful  if  there  is  drift  present  in  the  data.  One  can 
then  offset  the  data  by  the  average  of  the  first  and  last 
data  values.  This,  of  course,  works  only  if  both  points 
correspond  to  the  same  point  in  the  respiratory  cycle. 
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3)  Search  Che  data  and  determine  the  minimum  and  maximum  values 
of  the  data  for  each  channel.  Thus,  one  can  tell  at  a  glance 
the  excursions  of  the  chest  and  abdomen  and  also  the  volume 
excursion  for  that  data  set.  The  FORTRAN  listing  of  PLT 
is  in  Appendix  3. 


PROGRAM  COMPARE 

Program  COMPARE  prepares  a  data  set  for  plotting  one  channel 
data  versus  any  other  channel.  For  example,  one  can  plot  the  change  in 
anterior-posterior  (A-P)  chest  diameter  versus  volume,  or  one  can  compare 
the  change  in  A-P  chest  diameter  versus  the  change  in  lateral  chest 
diameter.  The  data  set  generated  is  compatible  with  the  automatic 
plotting  program  [1,100]  PLOT.  The  FORTRAN  listing  for  COMPARE  is  found 
in  Appendix  4. 

PROGRAMS  Ml  to  M7 

These  programs  determine  the  constants  for  the  seven  models  discussed 
below.  The  FORTRAN  listing  for  Model  7  is  found  in  Appendix  5.  For  the 
other  models,  the  array  "S"  is  altered  according  to  the  definitions 
in  reference  1. 

PROGRAM  CALC 

Program  CALC  generates  data  files  for  plotting  the  model  curves 
on  the  VERSATEC  plotter.  The  program  also  determines  whether  adding  a 
constant  to  the  predicted  volume  lowers  the  error.  The  program  asks 
the  user  whether  he  wishes  to  use  the  offset.  The  offset  should  be  used 
only  when  predictive  results  are  desired  and  not  for  calibration  curves. 
The  FORTRAN  listing  is  in  Appendix  6. 

The  equations  for  Models  1  through  7  are  given  below. 


i 
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Model  Is  V  -  Kj*(DC  -  DCO)  -  K2*(DA  -  DAO) 

Model  2:  V  -  K^CAPC  -  APCO)  -  K2*(APA  -  APAO) 

Model  3:  V  -  -  K2*(APA  -  APAO) 

Model  4:  V  -  Kj*(APC  -  APCO)  -  K2 

Model  5:  V  -  K^DC  -  K2* (DC/DA)  -  K3 

Model  6:  V  -  Kj*(DC  -  DCO)  -  K2*(DC/DA  -  DCO /DAO)  -  K3 


Model  7:  V  -  K^DC  -  K2*DA  -  K y 


where  V  is  volume,  DC  is  the  cross-sectional  area  of  the  chest,  DA  is  the 
cross-sectional  area  of  the  abdomen,  AFC  is  the  A-P  diameter  of  the 
chest,  and  APA  is  the  A-P  diameter  of  the  abdomen.  The  corresponding 

4* 


values  at  PRC  (measured  with  callipers)  are  followed  by  "0".  K, ,  K„ 


and  K3  are  unknown  constants  that  are  determined  by  minimising  the 
red  error. 
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In  the  previous  equations,  although  they  are  not  anatomical  models, 
we  still  assume  the  cross-sections  of  the  abdomen  and  chest  are  elliptical 
when  calculating  areas. 
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******************************* 


EXAMPLE  OF  USING  PROGRAM  ALLDAT  TO  INSPECT  DATA 


******************************************************************* 


RUN  ALLOAT 

FILE  NAME?  DKSPSRB.DAT I 1 
PRINT  EVERY  TTH  POINT?  10 
EVERY  10TH  POINT  IS  PRINTED 


0.5824 

0.2149 

0.4782 

-0.3183 

0.7462 

1.1353 

0.4739 

0.7131 

-0.0307 

0.8708 

1.4856 

0.5502 

0.9084 

0.2076 

0 . 8653 

1.3799 

0.5329 

0.9054 

0.1524 

0.7831 

1.0560 

0.3948 

0.7513 

-0.0481 

0.6531 

0.5427 

0.3056 

0.4279 

-0.2806 

0.5079 

0.0448 

0.0365 

0.2478 

-0.3968 

0.4805 

-0.2305 

-0.0843 

0 . 2066 

-0.5508 

0.3313 

0.1748 

0.0941 

0.2997 

-0.4259 

0.5216 

0.8534 

0.3948 

0.6079 

-0.1411 

0.6106 

1.3953 

0.5315 

0.8810 

0.1582 

0.7243 

1.4856 

0.5861 

0.9191 

0.2570 

0.7968 

1.1353 

0.4380 

0.7864 

0.0216 

0.6763 

0.6352 

0.30B5 

0.4858 

-0.2748 

0.5134 

0.1484 

0.1027 

0.2356 

-0.3823 

0.4395 

-0.1799 

-0.0814 

0.2021 

-0.4956 

0.3751 

-0.0235 

-0.0095 

0.2082 

-0.4840 

0.4381 

TT2  —  STOP 
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EXAMPLE  USING  PLT  TO  PREPARE  DATA  FOR  PLOTTING 

******************************************************************** 
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RUN  PLT 

INPUT  FILENAME  DK5PSRB.DAT f 1 
WRITE  MODIFIED  DATA  TO  FILETNO 
PTS  170  NO  OF  DATA  POINTS  FOR  DC  AYE  0.0 
0  0 

OFFSETS  0.000  0.000  0.000  0.000  0.000 

PREPARE  DATA  FOR  PLOT?  YES 

FILE  NAME  FOR  PLOT  DATA?  TEST .PLOf 1 

PLOT  EVERY  ?TH  DATA  POINT  2 


CHANNEL 

J 

MIN 

J 

MAX 

i.i  1 1  Kt! KhNl 

1 

155 

-0.250 

107 

1.538 

1 .78? 

2 

156 

-0.120 

106 

0.602 

0.722 

3 

159 

0.190 

104 

0.927 

0.737 

4 

155 

-0.562 

108 

0.266 

0 . 828 

5 

155 

0.298 

12 

0.883 

0.585 

ANY  MORE  ?  YES 
INPUT  FILENAME  DK JPSRB . DAT f 1 
WRITE  MODIFIED  DATA  TO  FILE?YES 
OUTPUT  FILE  NAME  PSRB.MODfl 
PTS  170  NO  OF  DATA  POINTS  FOR  DC  AVE  IrO 
1  0 


OFFSETS 

0.582 

0.215 

0.478 

-0.318 

0 . 746 

TRUNCATE  DATA?  NO 

PREPARE  DATA  FOR  PLOT?  YES 

FILE 

NAME 

FOR  PLOT 

DATA?  TEST 

•  PLOf  2 

PLOT 

EVERY 

?TH  DATA 

POINT  2 

CHANNEL  J 

MIN 

J 

MAX 

DIFFERFNC 

1 

155 

-0.833 

107 

0.956 

1.789 

2 

156 

-0.335 

106 

0.387 

0.722 

3 

159 

-0.288 

104 

0.449 

0.737 

4 

155 

-0.244 

108 

0.584 

0 . 828 

5 

155 

-0.448 

12 

0.137 

0.585 

ANY 

MORE 

?  NO 

TT2  —  STOP 
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****************************************************************** 


PREPARING  THE  PLOTS  FROM  DATA  SET  GENERATED  IN  PLT 


******************************************************************* 


GClflOOTPLOT 

>RUN  CllOfllOTPLOT.TSK 

PROGRAM  PLOT  09 ‘.05J 13  OB-AUG-79 

DATA  FILE?  C DEV .NAME . EXT » VERT  TEST.PL0J1 
NUMBER  OF  GRAPHS  ON  THIS  SHEET?  (MAX  10)  5 
X  MINIMUM  AND  X  MAXIMUM? 

Y  MINIMUM  AND  Y  MAXIMUM?  -l.»2. 

X  AXIS  LABEL?  <80  CHAR> 

SAMPLE  POINT (TIME ) 

Y  AXIS  LABEL?  <80  CHAR> 

OUTPUT ( L  *  CM ) 

NO.  OF  DECIMAL  DIGITS  FOR  LABELING  AXES?  2 
ANY  NOTES?  <YES  OR  N0>  YES 
NOTE?  <80  CHAR>  TEST. FLO »1 
POSITION?  (XfY)  ,5f 8. 

ANY  MORE?  NO 

- - CALCULATING  VECTORS - - 

WANT  GRID  PATTERN  ON  THE  GRAPH?  <YES  OR  N0>  NO 
091 08J 11  08-AUG-79 

ANY  MORE  PLOTS?  YES 

DATA  FILE?  CDEV . NAME . EXT i VERT  TEST.PL0»2 
NUMBER  OF  GRAPHS  ON  THIS  SHEET?  (MAX  10)  5 
X  MINIMUM  AND  X  MAXIMUM? 

Y  MINIMUM  AND  Y  MAXIMUM?  -l.»2. 

X  AXIS  LABEL?  <80  CHAR> 

SAMPLE  POINT ( TIME ) 

Y  AXIS  LABEL?  <80  CHAR> 

OUTPUT ( L » CM ) 

NO.  OF  DECIMAL  DIGITS  FOR  LABELING  AXES?  2 
ANY  NOTES?  <YES  OR  N0>  YES 
NOTE?  <80  CHAR>  TEST. FLO » 2 
POSITION?  ( X j Y )  .5*8. 

ANY  MORE?  YES 

NOTE?  <80  CHAR>  A  =  VOLUME*  B  *  AP  CHEST*  C  =  L.AT  CHEST 
POSITION?  ( X * Y )  .5*7.5 
ANY  MORE?  YES 

NOTE?  <80  CHAR>  D  =  AP  ABDOMEN  *  E  =  LAT  ABDOMEN 
POSITION?  ( X *  Y )  .5*7. 

ANY  MORE?  NO 

- - CALCULATING  VECTORS _ — - - 

WANT  GRID  PATTERN  ON  THE  GRAPH?  <YES  OR  N0>  NO 

09.*10:59  08-AUG-79 

ANY  MORE  PLOTS?  NO 

TT2  —  STOP  END  OF  PROGRAM  PLOT 

>*  WHAT  IS  YOUR  TERMINAL  NUMBER  *TTN‘  CST:  TT2 

>NFT  0PLTT2 

>PIP  PLTT2 . CMD *  *  *  PLTT2A . CMD i *  *  PLTT2B . CMD *  */DE 
>0  <EOF>  -8- 


******************************************************************* 


FITTING  VARIOUS  MODELS  TO  DATA  SET 


ft****************************************************************** 


RUN  Ml 


FILE  NAME  DN t PSRB ♦ DAT J 1 
INITIAL  VOLUME?  0 

APC*LATC,APAB,LATAB  =  ?  20.7,36. ,20.7,30.5 

APC  =  20.7  LC  =  36.0  APAB  *  20.7  LAB  =  30.5 

-0.24825090E+09 
V  =  A*<AC-ACO>  -  B*<AA  -  AAO ) 

A  =  0. 31541519E+02  B  =  -0.27743046E+02 

PRINT  EVERY  ?TH  PT?10 


POINT 

MEASURED 

CALC 

ERR 

1 

582.4 

559.3 

-23 .  1 

11 

1135.3 

1168.5 

33.2 

21 

1485.6 

1501.1 

15.4 

31 

1379.9 

1408.6 

28.7 

41 

1056.0 

1006.6 

-49.4 

51 

542.7 

534.7 

-8.0 

61 

44.8 

108.8 

64.0 

71 

-230.5 

-190.3 

40.3 

81 

174.8 

185.7 

10.9 

91 

853.4 

849.5 

-3.8 

101 

1395.3 

1371,7 

-23.6 

111 

1485.6 

1542.1 

56.4 

121 

1135.3 

1122.1 

-13.2 

131 

635.2 

573.8 

-61  .4 

141 

148.4 

153.5 

5.2 

151 

-179.9 

-133.6 

46.2 

161 

-23.5 

-30.4 

-7.0 

ERROR  SQ 

0.2928E+06 

ERRSGI/N 

0 . 1 723E+04 

CORRELATION 

COEFFICIENT 

0.9978 

ANY  MORE  ?N 
TT2  —  STOP 


0 • 4150E+02 
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RUN  M2 


FILE  NAME  DK J PSRB . HAT i 1 
INITIAL  VOLUME?  0 

ARC  t  LAT C  t  APAB » LATAB  =?  20 . 7 , 36 . , 20 . 7 » 30 . 5 

ARC  =  20.7  LC  =  36.0  ARAB  =  20.7  LAb  =  30.5 

-0 . 36518832E+03 

V  ~  A*<D1C-D1C0>  -  &*<L>1A  -  D1A0) 

f 

A  =  0 . 24866968E+04  B  =  -0 . 14099380E+03 

PRINT  EVERY  ?TH  RT?10 


POINT 

MEASURED 

CALC 

1 

582.4 

489.6 

11 

1135.3 

1174.2 

21 

1485.6 

1397.4 

31 

1379.9 

1346.7 

41 

1056.0 

974.9 

51 

542.7 

720 . 4 

61 

44.8 

34.9 

71 

-230.5 

-287.3 

81 

174.8 

173.9 

91 

853.4 

961.8 

101 

1395.3 

1343.9 

111 

1485.6 

1493.8 

121 

1135.3 

1092.1 

131 

635.2 

728 . 3 

141 

148.4 

201.6 

151 

-179.9 

-272.4 

161 

-23.5 

-91.8 

ERROR  SQ 

0. 7859E+06 

ERRSQ/N 

CORRELATION 

COEFFICIENT 

0.9935 

ANY  MORE  ?N0 
TT2  —  STOP 


ERK 

-92.7 
38.8 
-88 . 2 
-33.2 
-81 . 1 
177.6 
-9.9 
-56.8 

-0.9 

108.5 

-51.4 

8.2 

-43.2 

93.1 

53.2 
-92.5 
-68.4 

0 . 4623E+04  RMS  0.6799E+02 
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RUN  M3 


FILE  NAME  UK 5PSRB. BA T J 1 
INITIAL  VOLUME?  0 

ARC » LATCr ARAB fLATAB  =?  20 . 7 , 36 . , 20 . 7 , 30 . 5 

ARC  =  20 ♦ 7  LC  =  36.0  ARAB  =  20.7  LAB  =  30.5 

-0. 19503306E+04 
V  =  A*(l.)  -  B# (III A  -  D1A0) 

A  =  0. 10705438E+04  B  =  -0 . 22330249E+04 

PRINT  EVERY  ?TH  PT?10 


POINT 

MEASURED 

CALC 

ERR 

1 

582.4 

359 . 7 

-2 22.7 

11 

1135.3 

1002.1 

-133.3 

21 

1485.6 

1534.2 

48.6 

31 

1379.9 

1410.9 

31.0 

41 

1056.0 

963.1 

-92.9 

51 

542.7 

444.0 

-98.7 

61 

44.8 

184.5 

139.6 

71 

-230.5 

-159.4 

71.1 

81 

174.8 

119.6 

-55.2 

91 

853.4 

755.5 

-97.8 

101 

1395.3 

1423.9 

28.6 

111 

1485.6 

1644.5 

158.9 

121 

1135.3 

1118.9 

-16.5 

131 

635.2 

457.0 

-178.2 

141 

148.4 

216.9 

68.5 

151 

-179.9 

-36.2 

143.7 

161 

-23.5 

-10.2 

13.3 

ERROR  SQ 

0.2181E+07 

ERRSO/N 

0.1 283E+05 

correlation 

COEFFICIENT 

0.9815 

ANY  MORE  ?N0 
TT2  —  STOP 

> 


0. 1133E+03 


i 


.11. 


RUN  M4 


FILE  NAME  DK JPSRB . DAT i 1 
INITIAL  VOLUME?  0 

APC»LATC,APAB,LATAB  =?  20. 7, 36. ,20. 7,30. 5 

APC  =  20.7  LC  =  36.0  APAB  =  20.7  I  Ah  -  v,.r. 

-0. 14635099E+04 
V  *  AJ> ( D1C-D1C0 )  -  B*<1.) 

A  =  0. 26056 196E+04  B  =  0 . 51 137753E+02 

PRINT  EVERY  ?TH  PT?10 


POINT  MEASUREB 


1 

582.4 

11 

1135.3 

21 

1485.6 

31 

1379.9 

41 

1056.0 

51 

542.7 

61 

44.8 

71 

-230.5 

81 

174.8 

91 

853.4 

101 

1395.3 

111 

1485.6 

121 

1135.3 

131 

635,2 

141 

148.4 

151 

-179.9 

161 

-23,5 

CALC 

ERR 

508.9 

-73.4 

1183.7 

48.4 

1382.4 

-103.2 

1337.4 

-42.5 

977.5 

-78.5 

745.1 

202.4 

44.1 

-0.7 

-270.8 

-40 .  i 

194,0 

19.2 

977.5 

124.2 

1333.7 

-61.6 

1476.1 

-9.5 

1090.0 

-45.3 

752.6 

117.4 

216.5 

68.2 

-263.3 

-83.4 

-75.9 

-52 . 4 

ERROR  SQ  0 . 9401E+06  ERRSO/N  0.5530E+04  RMS 
CORRELATION  COEFFICIENT  0.9921 
ANY  MORE  ?N0 
TT2  —  STOP 

> 


0.74361+0 


-12- 


MPuPdata  FILE  NAME  UN J PSRB . DAT i 1 
APCrLATC»APAB»LATAB  =?  20 . 7 .  36 .  »  20 . 7  f  30 . 5 


APC  =  20.7 

LC  =  36.0 

ARAB  = 

20.7  l  Ah  = 

FRC  VOL  =  ? 

0 

FRC  VOLUME 

IN  L  «  0. 

000 

0 . 61679 439E+08  0 . 122921 75E+06  0 

. 10238236E+06 

0.1 22921 75E+06  0 . 24504753E+03  0 

. 204 101 35E+03 

0. 10238236E+06  0 . 20410135E+03  0 

. 1 7000000E+03 

0. 13622161E+05  0 . 75260568E+06  0 

•  742768.1 2E+08 

V  =  HO*AC  - 

VA*AC/AA  -  VC 

HO  0 . 55248625E+02  VA 

0 . 54526452E+04  VC 

PRINT  EVERY 

?TH  PTT10 

POINT 

MEASURED 

CALC 

ERR 

1 

582.4 

533 . 7 

-48  .  7 

11 

1 1 35 . 3 

1154.7 

19.3 

21 

1485.6 

1471.8 

- 1  3 . 8 

31 

1379.9 

1412.8 

32.9 

41 

1056.0 

1  030 .  1 

-  25 . 9 

51 

542 . 7 

573.9 

31 . 1 

61 

44.8 

60.3 

15.5 

71 

-230 . 5 

-202.3 

28.3 

81 

174.8 

169.5 

-5.3 

91 

853.4 

885 . 9 

32.5 

101 

1395.3 

1382.2 

-13.1 

111 

1485.6 

1527.5 

41.9 

121 

1135.3 

1 1  38 . 5 

3,1 

131 

635 . 2 

622 . 9 

-1  2.3 

141 

148.4 

130.8 

-17.5 

131 

-179.9 

-l 74 . 8 

5.  t 

161 

-23.5 

-63. 1 

-  >.9.6 

ERROR  SO 

0 . 1035E+06 

EKRSQ/N 

0 . O086E+03 

CORRELATION 

COEFFICIENT 

0.9991 

ANY  MORE  ?N0 
TT2  STOP 


M  )  .  5 


'o.71447678K+08 
0. 14035843E+06 
0.1 16S3972E+06 
0..  35471 72 7F+09 

0  ♦  2  A  0  3 9  7  7 0 E + 0  5 


RMS  0 . 2467E+0 


«*  -  r 


RUN  M6 

INPUT  DATA  FILE  NAME  DN : PSRB . DAT i 1 
APC»LATC> APAB.LATAB  =  ?  20. 7, 36 . . 20 . 7 , 30 . 5 

APC  =  20.7  LC  =  36.0  ARAB  =  20.7  LAD  = 

FRC  VOL  =  ?  0 

FRC  VOLUME  IN  L  =  0.000 

0.68712757E+05  0 . 60359670E+02  0 . 28849821E+04 
0.6 03596 70E +02  0 . 74073050E-01  0 . 34456037E+01 
0.28849821E+04  0 . 34456037E+01  0 . 1 7000000E+03 
0.13622166E+05  0 . 75260603E+06  0 . 74276729E+08 
-0.38364E-01  -0.28492E-04  -0.13270E-02 

V  =  H0*<AC-AC0>  -  VA# < AC/AA  -  ACO/AAO)  -  VC 
HO  0 . 5524863 1E+02  VA  0 . 54526372E+04  VC 


PRINT  EVERY 

?TH  PTT10 

POINT 

MEASURED 

CALC 

ERR 

1 

582.4 

533 . 7 

-48.7 

11 

1135.3 

1154.7 

19.3 

21 

1485.6 

1471.8 

-13.8 

31 

1379.9 

1412.8 

32.9 

41 

1056.0 

1030.1 

-25.9 

51 

542 . 7 

573.9 

31.1 

61 

44.8 

60.3 

15.5 

71 

-230.5 

-202.3 

28.3 

81 

174.8 

169.5 

-5.3 

91 

853.4 

885.9 

32.5 

101 

1395.3 

1382.2 

-13.1 

111 

1485.6 

1527.5 

41.9 

121 

1135.3 

1138.5 

3.1 

131 

635.2 

622.9 

-12.3 

141 

148.4 

130.8 

-17.5 

151 

-179.9 

-174.8 

5.1 

161 

-23.5 

-63.1 

-39.6 

ERROR  SO  0 

.1035E+06 

ERRSQ/N  0 

.  6086E+03 

NEW  OFFSET 

AND  VC  0. 

000  139.788 

CORRELATION 

COEFFICIENT 

0.9991 

REPRINT  DATA 

i  WITH  NEW  VC?  N 

ANY  MORE  ?N 
TT2  —  STOP 


30.5 


0 . 30638805E+07 
0 . 244924 16E+04 
0. 11683972E+06 
0. 19042154E+07 


0. 13978799E+03 


RMS  0 . 2467E+02 


■14' 


^  - 


RUN  H 7 

INPUT  DATA  FILE  NAME  DK J PSRB . DAT f 1 
APC»LATC>APAB»LATAB  =?  20 . 7 , 36 . , 20 . 7 , 30 . 5 

APC  *  20.7  LC  =  36.0  APAB  =  20.7  LAh  =  ■»>. 

PRC  VOL  =  ?  0 

FRC  VOLUME  IN  L  =  0.000 

0 . 616794 39E+08  0.51373341E+08  0 . 10238236E+06  0. 
0.5 137334 IE +08  0 . 42790040E+08  0 . 85276428E+05  0. 
0.10238236E+06  0 . 85276428E+05  0 . 1 7000000E+03  0. 
0 • 23749693E+10  0 . 1 05338 72E+ 12-0 . 31028363E+1 l  0. 
V  =  H0*AC  -  VA*AA  -  VC 

HO  0 • 443537 19E+ 02  VA  -0 . 1 3064743E +02  VC  0. 
PRINT  EVERY  ?TH  PT?10 


/  [  44  ' 

5947  7 
1  1683 
7737.,' 


(*1  '  HK  f-l/M 

5 t  7E  +  u8 
v  7  2 +'+06 

4vui-  +  l  4 


3  5  7 8 3 1  2 E + 0  3 


POINT 

MEASURED 

CAIC 

ERR 

1 

582.4 

533.7 

-48.  / 

11 

1135.3 

1154.7 

IV. 4 

21 

1485.6 

1472.3 

-13. 3 

31 

1379.9 

1412.5 

32.6 

41 

1056.0 

1029.8 

-26 . 2 

51 

542.7 

574.2 

3 1 . 5 

61 

44.8 

60.  1 

15.2 

71 

-230.5 

-201.9 

28 . 6 

81 

174.8 

169.6 

-5.2 

91 

853.4 

885.7 

33 . 4 

101 

1395.3 

1381.9 

-13.4 

111 

1485.6 

1527.7 

42.0 

121 

1135.3 

1138.1 

2.8 

131 

635 . 2 

623 . 3 

-12.0 

141 

148.4 

130.8 

-17.5 

151 

-179.9 

-175.1 

4.8 

161 

-23.5 

-63.2 

-39. 7 

tROR  SQ 

0 ♦ 1034E+06 

ERRSQ/N 

0  •  6080E+03 

CORRELATION  COEFFICIENT 
ANY  MORE  ?N0 
TT2  —  STOP 


0.9991 


0 . 2466+  +0 


******************************************************************* 


PREPARING  PLOTS  FOR  THE  CURVE  FITS 


*****************************************************  *************** 


RUN  CALC 

DO  YOU  NEED  MODEL  DESCRIPTION?  YES 


MODEL  CALCULATIONS  FOR  ELIPTIC  CYLINDERS 


TYPE  1  V  =  A*  <  AC-ACO )  -  B# ( AAB-AABO ) 

TYPE  2  V  =  A*  ( D1C  <  DI  AM)  -D1CO )  -  B*  <  III  A  (  DI  AM  >  -  lMAO) 

TYPE  3  V  =  A  -  B#(APAB(DIAM)-  APABO) 

TYPE  4  V  *  A*  ( AF'C  ( DI  AM )  -AF'CO )  -  B 

TYPE  5  V  =  A*  <  AC )  -  B*  <  AC/AAB )  -  C 

TYPE  6  V  ~  A* (AC-ACO)  -  B* < *AC/AAB-ACO/AARO >  -  C 

TYPE  7  V  »  A*AC  -  B*AAB  -  C 


INPUT  DATA  FILE  NAME  DK  J F'SRB .  DAT !  1 
MODEL  TYPE?  I 


COMPARISON  FOR  MODEL  TYPE  1 


PLOT  COMPARISON  ?  YES 
FILE  NAME?  M1.PLOI1 

APC»LATC»APAB»LATAB  =?  20.7.36. >20.7,30.5 

APC  =  20.7  LC  =  36.0  APAB  =  20.7  LAB  =  30.5 

FRC  VOL  =  ?  0 

FRC  VOLUME  IN  L  =  0.000 

A>B>C  =?  31 .541 51 9 >-27. 743046 

A  0.31 54151 9E+02  B  -0.27743046E+02  C  O.OOOOOOOOF+OQ 
PRINT  EVERY  ?TH  PTT169 


POINT 

MEASURED 

CALC 

ERR 

1 

582.4 

559.3 

-23.1 

170 

639.7 

697.9 

58.2 

ERROR  SO 
OFFSET  ? 

0.2928E+06 

NO 

ERRSQ/N 

0. 1722E+04  RMS 

0 » 4150E+02 

CORRELATION  COEFFICIENT 

0.9978 

-16  - 

REPRINT  DATA  WITH  NEW  C?  NO 
EVERY  TTH  POINT  PLOTTED?  2 


ANY  MORE  ON  1HIS  FILE?  YES 
MODEL  TYPE?  2 


COMPARISON  FOR  MODEL  TYPE  2 


PLOT  COMPARISON  ?  YES 
FILE  NAME?  M2.PL0.1 
NEW  CHEST  DIMENSIONS  ?NO 

FRC  VOL  =  ?  0  ! 

FRC  VOLUME  IN  L  =  0.000 

A » B » C  =?  2486.6968.-141 .99380 

A  0 . 24866968E+04  B  -0 . 1 4 1 99380E+03  C  0 . OOOOOOOOE+OO 
PRINT  EVERY  ?TH  PT?169 


POINT 

MEASURED 

CALC 

ERR 

1 

582.4 

489.3 

-93.  t 

170 

639 . 7 

677 . 3 

37.6 

ERROR  SO 

nrrprr  o 

0 . 7859E+0A 

un 

ERRSG/N 

0 • 4623E+04  RMS 

0 • 6799E+02 

OFFSET  ?  NO 

CORRELATION  COEFFICIENT  0.9935 
REPRINT  DATA  WITH  NEW  C?  NO 
EVERY  ?TH  POINT  PLOTTED?  2 
ANY  MORE  ON  THIS  FILE?  YES 
MODEL  TYPE?  3 


COMPARISON  FOR  MODEL  TYPE  3 


PLOT  COMPARISON  ?  YES 

FILE  NAME?  M3.PL051 

NEW  CHEST  DIMENSIONS  ?N0 

FRC  VOL  =  ?  N\N\0 

FRC  VOLUME  IN  L  =  0.000 

A.B.C  =?  1070 . 5438  . -2233 . 0249 

A  0. 10705438E+04  B  -0 . 22330249E+04  C  0 . OOOOOOOOE+OO 
PRINT  EVERY  ?TH  PT?169 


POINT 

MEASURED 

CALC 

ERR 

1 

582.4 

359.7 

-222.7 

170 

639.7 

671 . 1 

31.5 

ERROR  SQ 
OFFSET  ? 

0.2181E+07 

NO 

ERRSG/N 

0. 1283E+05  RMS 

0 . 1 1 33F+03 

CORRELATION  COEFFICIENT  0.9815 
REPRINT  DATA  WITH  NEW  C?  NO 
EVERY  ?TH  POINT  PLOTTED?  2 
ANY  MORE  ON  THIS  FILE?  YES 
MODEL  TYPE?  4 


COMPARISON  FOR  MODEL  TYPE  4 


PLOT  COMPARISON  ?  Y 
FILE  NAME?  M4.PL0I1 
NEW  CHEST  DIMENSIONS  ?N 
FRC  VOL  »  ?  0 

FRC  VOLUME  IN  L  *  0.000 

A»BfC  -?  2605.619651.137753 

TT7  —  FXTTTNA  nilF  TO  FRROR  64 


ffrRMkED  MODEL  DESCRIPTION?  NO 
INPUT  DATA  FILE  NAME  DKJPSRB.DATJ 1 
MODEL  TYPE?  4 


COMPARISON  FOR  MODEL  TYPE  4 


PLOT  COMPARISON  ?  YES 
FILE  NAME?  M4.PLOU 

APC.LATC.APAB.LATAB  '  =?  20 . 7, 36 . .20 . 7 , 30 . 5 

APC  *  20.7  LC  =  36.0  APAB  =  20.7  LAB  *  30.5 

FRC  VOL  *  ?  0 

FRC  VOLUME  IN  L  =  0.000 

A.B.C  =?  2605.61 96 .51 • 137753 

A-  0 .26056196E+04  B  0.51137753E+02  C  0 .OOOOOOm/t-  +00 

■PRINT  EVERY  ?TH  PTT169 


.POINT  MEASURED  CALC  ERR 

1  582.4  508.9  -73.4 

170  639.7  685.1  45.5 

ERROR  SO  0.9401E+06  ERRSQ/N  0.5530E+04  RMS  0.74.i>6h+07 

OFFSET  ?  NO 

CORRELATION  COEFFICIENT  0.9921 
REPRINT  DATA  WITH  NEW  C?  NO 
EVERY  ?TH  POINT  PLOTTED?  2 
ANY  MORE  ON  THIS  FILE?  YES 
MODEL  TYPE?  5 


COMPARISON  FOR  MODEL  TYPE  5 


PLOT  COMPARISON  ?  Y 

FILE  NAME?  M5.PL0J1 

NEU  CHEST  DIMENSIONS  ?N0 

FRC  VOL  *  ?  0 

FRC  VOLUME  IN  L  =  0.000 

A.B.C  -?  55.248625.5452,6452.26039.72 

A  0. 55248625E+02  B  0.54526452E+04  C  0.26039720E+05 
PRINT  EVERY  ?TH  PT?169 


POINT  MEASURED  CALC  ERR 

1  582.4  533.7  -48.7 

170  639.7  634.3  -5.4 

ERROR  SO  0. 1035E+06  ERRSO/N  0.6086E+03  RMS  0.2467E+0, 

OFFSET  ?  NO 

CORRELATION  COEFFICIENT  0.9991 
REPRINT  DATA  WITH  NEU  C?  NO 
EVERY  ?TH  POINT  PLOTTED?  2 
ANY  MORE  ON  THIS  FILE?  Y 
MODEL  TYPE?  6 


COMPARISON  FOR  MODEL  TYPE  6 


PLOT. COMPARISON..!.-  Y _ 


/  - 


FILE  NAME?  M6.PL0K 

NEW  CHEST  DIMENSIONS  ?N 

FRC  VOL  =  ?  0 

FRC  VOLUME  IN  L  =  0.000 

A,B,C  =?  55 . 248631 .5452.6372  *  139 . 78799 

A  0. 5524863 1E+02  B  0.54526372E+04  C 

PRINT  EVERY  ?TH  PT?169 


0.13978799E+03 


? 

t. 

i 


POINT  MEASURED  CALC 

1  582.4  533.7 

170  639.7  634.3 


ERR 

-48.7 

-5.4 


ERROR  SO  0 . 1035E+06  ERRSO/N  0.6086E+03  RMS  0.2467E+02 
OFFSET  ?  NO 

CORRELATION  COEFFICIENT  0.9991 
REPRINT  DATA  WITH  NEW  C?  NO 
EVERY  ?TH  POINT  PLOTTED?  2 
ANY  MORE  ON  THIS  FILE?  YES 
MODEL  TYPE?  7 


COMPARISON  FOR  MODEL  TYPE  7 


PLOT  COMPARISON  ?  YES 

FILE  NAME?  M7.PL0, 1 

NEW  CHEST  DIMENSIONS  ?N0 

FRC  VOL  *  ?  0 

FRC  VOLUME  IN  L  =  0.000 

A.B.C  =?  44.353719,-13.064743,32578.312 

A  0.44353719E+02  B  -0 . 13064743E+02  C  0.32578312E+05 
PRINT  EVERY  ?TH  PT?169 


POINT  MEASURED  CALC  ERR 

1  582,4  533.7  -48.7 

170  639.7  634.6  -5.1 


ERROR  SQ  0, 1034E+06  ERRSQ/N  0.6080E+03  RMS  0.2466E+02 
OFFSET  ?  NO 

CORRELATION  COEFFICIENT  0.9991 
REPRINT  DATA  WITH  NEW  C?  NO 
EVERY  ?TH  POINT  PLOTTED?  2 
ANY  MORE  ON  THIS  FILE?  NO 
ANY  MORE  ?N0 
TT2  —  STOP 
> 


/  •  - 


****************************************************************** 


PLOTTING  THE  CURVE  FITS 


******************************************************************* 


(?C  1 » 1003PL0T 

>RUN  CllO. 1103PL0T.TSK 

PROGRAM  PLOT  10S01M1  08-AUG-79 

DATA  FILE?  C DEV. NAME . EXT » VERT  Ml.FLO.l 
NUMBER  OF  GRAPHS  ON  THIS  SHEET?  (MAX  10)  2 
X  MINIMUM  AND  X  MAXIMUM? 

Y  MINIMUM  AND  Y  MAXIMUM? 

X  AXIS  LABEL?  <80  CHAR> 

SAMPLE (TIME 

Y  AXIS  LABEL?  <80  CHAR  > 

VOLUME 

NO.  OF  DECIMAL  DIGITS  FOR  LABELING  AXES?  0 
ANY  NOTES?  <YE3  OR  N0>  YES 
NOTE0  <80  CHAR>  M1.PL0J1 
POSITION?  ( X » Y  )  ,5»8, 

ANY  MORE?  NO 

- - CALCULATING  VECTORS - - 

WANT  GRID  PATTERN  ON  THE  GRAPH?  <YES  OR  N0>  NO 
10:02J42  08-AUG-79 

ANY  MORE  PLOTS?  YES 

DATA  FILE?  CDEV . NAME . EXT  5  VERT  M2.PL0J1 
NUMBER  OF  GRAPHS  ON  THIS  SHEET?  (MAX  10)  2 
X  MINIMUM  AND  X  MAXIMUM? 

Y  MINIMUM  AND  Y  MAXIMUM? 

X  AXIS  LABEL?  <80  CHAR> 

SAMPLE (TIME) 

Y  AXIS  LABEL?  <80  CHAR> 

VPO\QP\OLUME 

NO.  OF  DECIMAL  DIGITS  FOR  LABELING  AXES?  0 
ANY  NOTES?  <YES  OR  N0>  YES 
NOTE?  <80  CHAR>  M2.PL0J1 
POSITION?  (X»Y)  .5»8. 

ANY  MORE?  NO 

- __ ..CALCULATING  VECTORS___ - 

WANT  GRID  PATTERN  ON  THE  GRAPH?  <YES  OR  N0>  NO 
10:03M0  08-AUG-79 

ANY  MORE  PLOTS?  YES 

DATA  FILE?  C DEV .NAME. EXT *VER3  M3.PL0J1 
NUMBER  OF  GRAPHS  ON  THIS  SHEET?  (MAX  10)  2 
X  MINIMUM  AND  X  MAXIMUM? 

Y  MINIMUM  AND  Y  MAXIMUM? 

X  AXIS  LABEL?  <80  CHAR> 

SAMPLE (TIME) 

Y  AXIS  LABEL?  <80  CHAR> 

VOLUME 

NO.  OF  DECIMAL  DIGITS  FOR  LABELING  AXES?  0 
ANY  NOTES?  <YES  OR  N0>  YES  _20_ 

NOTF*  <R0  CHAR>  M3. PI  0*1 


POSITION?  (  X  .  Y  )  .5.8. 

ANY  MORE7  NO 

- - CALCULATING  VECTORS - - 

WANT  GRID  PATTERN  ON  THE  GRAPH?  <YES  OR  NO  NO 
10:04J30  08-AUG-79 

ANY  MORE  PLOTS7  YES 

DATA  FILE?  t DEV . NAME . EXT » VERI  M4.PL0.1 
NUMBER  OF  GRAPHS  ON  THIS  SHEET?  (MAX  10)  2 
X  MINIMUM  AND  X  MAXIMUM? 

Y  MINIMUM  AND  Y  MAXIMUM7 
X  AXIS  LABEL?  <80  CHAR> 

SAMPLE (TIME) 

Y  AXIS  LABEL?  <80  CHAR,* 

VOLUME 

NO.  OF  DECIMAL  DIGITS  FOR  LABELING  AXES7  0 
ANY  NOTES?  <YES  OR  N0>  YES 
NOTE?  <80  CHAR>  M4.PL0.1 
POSITION?  ( X » Y )  .5.8. 

ANY  MORE?  NO 

- - CALCULATING  VECTORS - - 

WANT  GRID  PATTERN  ON  THE  GRAPH7  <YES  OR  NO.-  NO 
10:05:24  08-AUG-79 

ANY  MORE  PLOTS?  YES 

DATA  FILE?  CDEV . NAME . EXT i VERT  M5.PL0J1 
NUMBER  OF  GRAPHS  ON  THIS  SHEET7  (MAX  10)  2 
X  MINIMUM  AND  X  MAXIMUM? 

Y  MINIMUM  AND  Y  MAXIMUM7 
X  AXIS  LABEL?  <80  CHAR> 

TI\IT\SAMF"LE  <  TIME ) 

Y  AXIS  LABEL?  <80  CHAR> 

VOLUME 

NO.  OF  DECIMAL  DIGITS  FOR  LABELING  AXES?  0 
ANY  NOTES?  <YES  OR  N0>  YES 
NOTE?  <80  CHAR>  M5.PL0.1 
POSITION?  ( X » Y )  .5.8, 

ANY  MORE?  NO 

- - CALCULATING  VECTORS - - 

WANT  GRID  PATTERN  ON  THE  GRAPH?  <YES  OR  N0>  NO 
10:06:20  08-AUG-79 

ANY  MORE  PLOTS?  YES 

DATA  FILE?  CDEV .NAME . EXT  5  VERT  M6.PL0J1 
NUMBER  OF  GRAPHS  ON  THIS  SHEET?  (MAX  10)  2 
X  MINIMUM  AND  X  MAXIMUM? 

Y  MINIMUM  AND  Y  MAXIMUM? 

X  AXIS  LABEL?  <80  CHAR> 

SAMPLE (TIME) 

Y  AXIS  LABEL?  <80  CHAR> 

VOLUME 

NO.  OF  DECIMAL  DIGITS  FOR  LABELING  AXES?  0 
ANY  NOTES?  <YES  OR  N0>  YES 
NOTE?  <80  CHAR>  M6.PL0J1 
POSITION?  (X»Y)  .5.8. 

ANY  MORE?  NO 

- - CALCULATING  VECTORS - - 

WANT  GRID  PATTERN  ON  THE  GRAPH?  <YES  OR  NO,-  NO 
10:07U0  08-AUG-79 

ANY  MORE  PLOTS?  YES 

DATA  FILE?  CDEV. NAME. EXT JVER3  M7.PL0J1 
NUMBER  OF  GRAPHS  ON  THIS  SHEET?  (MAX  10)  7 
X  MINIMUM  AND  X  MAXIMUM? 

Y  MINIMUM  AND  Y  MAXIMUM? 

X  AXIS  LABEL?  <80  CHAR> 

SAMPLE (TIME)  -21- 

Y  AXIS  LABEL?  <80  CHAR> 

I  VOLUME 

wn.  OF  nFr.TMAt  titottr  FOR  I  ARFl  TNG  AXFR7  0 


ANY  NOTES?  <YES  OR  NO>  YES 
NOTE?  <80  CHAR>  M7.PL051 
POSITION?  ( X » Y )  .5*8. 

ANY  MORE?  NO 

- - CALCULATING  VECTORS - - 

WANT  GRID  PATTERN  ON  THE  GRAPH7  <YES  OR  N0>  NO 
10507:58  08-AUG-79 

ANY  MORE  PLOTS?  NO 

TT2  —  STOP  END  OF  PROGRAM  PLOT 

>*  WHAT  IS  YOUR  TERMINAL  NUL3ER  *  TTN  *  CS3  5  TT2 

>NFT  0PLTT2 

>PIP  PLTT2. CMD5  PLTT2A . CMD i *  » PLTT2B . CMD *  */DE 
>e  <EOF> 

> 
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C0RE-19K.  UIOC1AO»13  ALLDAT »ALLDAT/-SP=ALLDAT 

C*******************M*************M********»********************* 

C 

c 

C  ALLDAT 

C 

C 

C **************************************************************** 

C  THIS  PROGRAN  PRINTS  THE  RAW  DATA  THAT  HAS  BEEN  SAMPLED  WITH 

C  DATLOG.  IT  HAY  BE  USED  TO  GET  A  FEELING  FOR 

C  THE  DATA  BUT  PRINTING  EVERY  PIECE  OF  DATA  WASTES  MUCH  TIME. 

C  YOU  MAY  ASK  FOR  ONLY  REPRESENTATIVE  POINTS  TO  BE  PRINTED 

C******************************************************************** 


0001 

DIMENSION  A<6) »B<200) . FILE ( 10) 

0002 

30 

FORMAT  <5<3X»F9.4)) 

0003 

700 

TYPE  29 

0004 

29 

FORMAT < '  FILE  NAME?  ',♦> 

0003 

28 

FORMAT (6A4) 

0006 

ACCEPT  55 .FILE 

0007 

33 

FORMAT <10A4) 

0008 

CALL  ASSIGN <3 .FILE. 0) 

0009 

CALL  FDBSET ( 3 . ' READONLY ' ) 

0010 

READ ( 3 . ERR™700 ) DUM 

0011 

REWIND  3 

0012 

TYPE  705 

0013 

705 

FORMAT < '  PRINT  EVERY  ?TH  POINT?  '.»> 

0014 

ACCEPT  706. N 

0015 

706 

FORMAT (13) 

0016 

TYPE  707. N 

0017 

707 

FORMAT < '  EVERY  ',I3.'TH  POINT  IS  PRINTED'///) 

0018 

IC  -  N-l 

0019 

1 

CONTINUE 

0020 

READ  (3.END-50)  < A< I >  » I“1 .6 ) 

0021 

IC  -  IC  +  1 

0022 

IF< IC/N4N.NE. IC)GO  TO  1 

0024 

TYPE  30»<A(K)»K=1.5) 

0023 

00  TO  1 

0026 

98 

CONTINUE 

0027 

31 

FORMAT <2X» 13) 

0028 

30 

STOP 

0029 

END 
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0001 

0002 

0003 

0004 

0005 

0000 

0007 

0008 

0009 

0010 

0012 

0013 

0014 

0015 

0010 

0017 

0018 

0019 

0020 

0021 

0022 

0023 

0024 

0020 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0030 

0037 


C************************************************************* 

c 

c 

C  PLT 

C 

C 

c********************»************************«***************** 

C  THIS  IS  THE  MAIN  DATA  REDUCTION  PROGRAM  FOR  DATA  FROM 

C  DATLOG.  THE  PROGRAM  WILL  TELL  YOU  MINIMUM  AND 

C  MAXIMUM  VALUES  FOR  EACH  CHANNEL , AND  PREPARE  THE  DATA  FOR 

C  eClrlOO]  PLOT.  YOU  MAY  ALSO  'DC  OFFSET*  THE  DATA 

C  BY  THE  AVERAGE  OF  THE  FIRST  *N1*  AND  LAST  *N2* 

C  DATA  POINTS.  IF  YOU  WANT  THE  BEGINING  DO  BE  ZERO  USE  N1>0 

C  AND  N2-0.  IF  YOU  WANT  THE  LAST  TO  BE  ZERO  USE  N1=0  AND  N2>0 

C  YOU  MAY  OF  COURSE  USE  N1  -  N2  *  0  AND  YOU  WILL  GET  NO 

C  OFFSET.  THE  DATA  AFTER  BEING  OFFSET  WILL  BE  WRITTEN  TO  A 

C  NEW  FILE  IF  YOU  CHOOSE.  REMEMBER  YOU  MAY  ONLY  PLOT  100 

C  DATA  POINTS  IN0C1, 1003PL0T. 

C*»***********************«*M****************»**********«****** 

DIMENSION  SUM<50),JMIN(5>,JMAX<5> 

DIMENSION  X<A),XMIN<5),XMAX<5),Y<5,1250) 

DIMENSION  FILE< 10) 

DATA  YB/1HY/ 

GO  TO  700 
751  CONTINUE 
TYPE  755 

755  FORMAT <'  ANY  MORE  T  ',•> 

ACCEPT  721 »RRR 
IF  <RRR.NE.YB) GO  TO  752 

700  TYPE  701 

701  FORMAT <'  INPUT  FILENAME  ',*> 

ACCEPT  702, FILE 

702  FORMAT <10A4) 

CALL  A8SIGN(3,FILE,0) 

CALL  FDBSET <  3 , ' READONLY ' ) 

READ  <  3 , ERR" 7 00 ) DUM 
REWIND  3 
TYPE  720 

720  FORMAT < '  WRITE  MODIFIED  DATA  TO  FILE?' »•> 

ACCEPT  721, R 

721  FORMAT (Al) 

IF(R.NE.YB)  00  TO  725 

710  TYPE  711 

711  FORMAT < '  OUTPUT  FILE  NAME  ',♦) 

712  ACCEPT  702, FILE 
CALL  A88IGN<4,FILE,0> 

CALL  FDB8ET  <  4 , ' NEW ' ) 

725  IC  -  0 

1  READ (3, END-50)  <X<I)»I-1,6> 

IC  -  IC  ,+  1 
DO  2  J  •  1,5 

2  Y( J,IC)  -  X<J> 

00  TO  1 

50  TYPE  51, IC 
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0038 

51 

FORHAT < '  PTS  '»I4»'  NO  OF  DATA  POINTS 

0039 

ACCEPT  52 f  N1 *N2 

0040 

52 

FORMAT! 13. 13) 

0041 

TYPE  296 »N1 »N2 

0042 

296 

FORMAT  !  3X  f  2 1 3 ) 

0043 

DO  4  J  -  If 5 

0044 

SUM ( J )  -  0 

0045 

IF(Nl.LE.O)  GO  TO  6 

0047 

DO  5  K  -  lfNl 

0048 

SUM! J >  -  SUM ! J )  +Y ! J r K > 

0049 

5 

CONTINUE 

0050 

6 

IF!N2.LE.O)  GO  TO  8 

0052 

DO  7  K  >  1»N2 

0053 

SUM! J)  -  SUM! J)  +  Y! Jf IC-K+1) 

0054 

7 

CONTINUE 

0055 

8 

CONTINUE 

0056 

IF!N1+N2.NE.0)SUM!J)  »  SUM ! J ) / ! N 1 +N2 ) 

0058 

DO  10  K  -  lfIC 

0059 

10 

Y! J»K)  -  Y! JrK)  -  SUM!J) 

0060 

4 

CONTINUE 

0061 

TYPE  58f !SUM!J)» J-1.5) 

0062 

58 

FORMAT!'  OFFSETS  ' t 5!2X»F8.3> ) 

0063 

IF!R.NE.YB>GO  TO  731 

0065 

NSTA  -  1 

0066 

NSTO  -  IC 

0067 

TYPE  1300 

0068 

1300 

FORMAT!'  TRUNCATE  DATA?  'r*> 

0069 

ACCEPT  721 fR4 

0070 

IF!R4.EQ. YB>  TYPE  1310 

0072 

IF!R4.EQ.YB)  ACCEPT  1305 f NSTA f NSTO 

0074 

1310 

FORMAT!'  1ST  AND  LAST  POINTS  T  'f»> 

0075 

1305 

FORMAT! 214) 

0076 

DO  730  K«  NSTA f NSTO 

0077 

URITE!4)!Y!JfK)f J-lf5)fDUM 

0078 

730 

CONTINUE 

0079 

731 

CONTINUE 

0080 

TYPE  740 

0081 

740 

FORMAT!'  PREPARE  DATA  FOR  PLOT?  'r*> 

0082 

ACCEPT  721 fRR 

0083 

IF!RR.NE.YB)GO  TO  750 

0085 

TYPE  760 

0086 

760 

FORMAT!'  FILE  NAME  FOR  PLOT  DATAT  'r*> 

0087 

ACCEPT  702 » FILE 

0088 

CALL  A88I8N! 1 f FILE f 0 ) 

0089 

CALL  FDBSETUf  'NEW') 

0090 

TYPE  60 

0091 

60 

FORMAT!'  PLOT  EVERY  TTH  DATA  POINT  'f 

0092 

77 

FORMAT !2XrFS. Of ' t ' fF9.4) 

0093 

72 

FORMAT! I3f'f 'fll»'»') 

0094 

ACCEPT  61 f N3 

0095 

61 

FORMAT! 12) 

0096 

NC-  ZC/N3 

0097 

SF!NC«N3fME.ZC>  NC-NC  4  1 

0099 

DO  70  1  -  IrS 

-27- 
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0100 

0101 

0102 

0103 

0104 

75 

0105 

70 

0106 

750 

0107 

0108 

0109 

0110 

0111 

0112 

0113 

0115 

0117 

0119 

0121 

90 

0122 

0123 

800 

0124 

0125 

89 

0126 

0127 

0129 

95 

0130 

0132 

0133 

752 

0134 

WRITE ( 1 p72>  NC»  I 
DO  73  K  -  1 » IC.N3 
A-K 

WRITE! 1*77)  A.Y(I.K) 

CONTINUE 
CONTINUE 
CONTINUE 
DO  90  J-1.5 
XMIN(J)  ■  Y< J* 1 ) 

XHAX  !J>“Y!J»1> 

JMIN(J)  -  1 
JMAX(J)  =  1 
DO  90  K  «  l.IC 

IF!XHIN(J> .GT.Y! J.K) >  JMIN(J>  -  K 
IF ( XHIN( J > .GT.Y! J.K) >XHIN!J>«  Y(J*K> 

IF (XHAX( J) .LT.Y!J»K) ) JHAX( J)  -  K 
IF ( XHAX !  J  )  ♦  LT  •  Y  !  J  »  K  )  >  XMAX <J)*Y!J»K> 

CONTINUE 
TYPE  800 

FORHAT  !  4X  »  '  CHANNEL  J  HIN  J  HAX 

1 DIFFERENCE  ' ) 

DO  89  K  -  1.5 

TYPE  95.K.JHIN!K).XHIN!K>.JHAX!K>.XHAX!K).XHAX(K)-XHIN<K> 
CALL  CLOSE ( 3 ) 

IF ! RR . EG  »  YB ) CALL  CLOSE! 1) 

FORHAT !7X.I3.4X.I4. 3X  »  F8 • 3  »  3X  » 1 4  »  3X  »  F8 • 3  »  5X . FB . 3  > 

IFCR.EQ. YB)CALL  CLOSE! 4 > 

GO  TO  751 

STOP 

END 
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C********************************************«***************** 

c 

c 

C  COMPARE 

C 

C 

C************************************************************** 

C 

C  THIS  PROGRAM  COMPARES  THE  OUTPUT  OF  ANY  TWO  CHANNELS  OF 

C  DATA  IN  A  FORMAT  CONSISTENT  WITH  THE  OUTPUT  OF  DATLOG 

C 
C 

C  **M*»**********«*************************«******************** 


0001 

DIMENSION  B<6> 

0002 

DIMENSION  X<6.999> 

0003 

DIMENSION  FILE(IO) 

0004 

DATA  Y/1HY/ 

0005 

100 

FORMAT  < 10A4  > 

0006 

110 

FORMAT !13»'»'»Il?'i') 

0007 

120 

FORMAT <2XrF8. 4 » ' * '  »F8.4> 

0008 

135 

CONTINUE 

0009 

NN  -  0 

0010 

1 

TYPE  101 

0011 

101 

FORMAT !  // '  INPUT  FILE?  '»*> 

0012 

ACCEPT  IOOpFILE 

0013 

CALL  AS8IGN<3»FILE»0) 

0014 

CALL  FDB8ET<3f 'READONLY' > 

0015 

READ ( 3  f  ERR* 1 ) DUM 

0016 

REWIND  3 

0017 

TYPE  102 

0018 

102 

FORMAT < '  OUTPUT  FILET  ' »•) 

0019 

ACCEPT  100'FILE 

0020 

CALL  A8SIGN!4»FILE*0> 

0021 

CALL  FDB8ET!4r'NEW'> 

0022 

TYPE  103 

0023 

103 

FORMAT < '  CHANNEL  T  -  X?  ' »*> 

0024 

ACCEPT  104vNl 

0025 

104 

FORMAT  < 1 1 ) 

0026 

TYPE  105 

0027 

105 

FORMAT! '  CHANNEL  ?  -  YT  '»♦> 

0028 

ACCEPT  104 tN2 

0029 

IC  -  0 

0030 

READ! 3) (B(I)»I«1»6) 

0031 

XMIN  -  B!N1 ) 

0032 

XMAX  -  XMIN 

0033 

YMIN  -  B!N2) 

0034 

YMAX  -  YMIN 

0035 

REWIND  3 

0034 

25 

READ!3»END-50>!X!I'IC+1)»I-1'6> 

•037 

IF!X!N1 r IC+1 > »LT •XMIN)  XMIN  -  X!N1»IC+1> 

0039 

IF!X!N1»IC91>*GT«XMAX)  XMAX  -  X!N1'IC+1) 

0041 

IF(X<N2»IC+1).LT.YMIN>  YMIN  -  X!N2»IC+1> 

0043 

IF(X!N2f 1C+1) t8T> YMAX)  YMAX  »  X!N2tIC+l> 

0045 

XC  -  ZC  ♦  1 

-30- 
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0046 

GO  TO  25 

0047 

50 

WRITE (4rllO)ICf  NN 

0048 

TYPE  150'XMIN'XMAX 

0049 

TYPE  151 »  YHIN.YMAX 

0050 

150 

FORMATS  XHIN  '»F8.4*'  XHAX 

'  »F8»4) 

00S1 

151 

FORMAT <'  YMIN  ' »F8. 4» '  YMAX 

'  »F8.4> 

0052 

DO  301  -lfIC 

0053 

30 

URITE(4»120>X(Nlf I)fX(N2rI) 

0054 

CALL  CLOSE< 3) 

0055 

CALL  CLOSE ( 4  > 

0056 

TYPE  130 

0057 

130 

format:'  any  more?  ',*> 

0058 

ACCEPT  131 fYB 

0059 

if:yb.eq.y>  go  to  135 

0061 

131 

FORMAT :A1> 

0062 

STOP 

0063 

END 

FORTRAN  IV  STORAGE  HAP 


NAME 

OFFSET 

ATTRIBUTES 

B 

000006 

REALM 

ARRAY  :6> 

X 

000036 

REALM 

ARRAY  :6»999>  VECTORED 

FILE 

056706 

REALM 

ARRAY  :iO) 

Y 

056756 

REALM 

VARIABLE 

NN 

063222 

INTEGER* 2 

VARIABLE 

A8SIGN 

000000 

REALM 

PROCEDURE 

FDB8ET 

000000 

REALM 

PROCEDURE 

DUN 

063224 

REALM 

VARIABLE 

N1 

063230 

INTEGERB2 

VARIABLE 

N2 

063232 

INTEGERB2 

VARIABLE 

IC 

063234 

INTEGER* 2 

VARIABLE 

I 

063236 

INTEGER*2 

VARIABLE 

XM1N 

063240 

REALM 

VARIABLE 

XHAX 

063244 

REALM 

VARIABLE 

YHIN 

063250 

REALM 

VARIABLE 

YMAX 

063254 

REALM 

VARIABLE 

CL08E 

000000 

REALM 

PROCEDURE 

YB 

063260 

REALM 

VARIABLE 
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C****************************************************************** 

c 

c 

C  M5  CURVE  FIT  FOR  MODEL  5  3  PARAMETER  MODEL 

C 

C 

C*****************************M*********************************** 

0001  DOUBLE  PRECISION  DC  *  DA  r  V  r  APC  ? XLC » APAB  » XLAB  » VFRC 

0002  DOUBLE  PRECISION  PI»D1C»D2C'D1A'D2A'H0'VA'VC'SER'SX 

1 »SX2»SYrSY2»SXY  > VM»D2D» VV  »ER»CC 
0003  DIMENSION  FILE (10) 

0004  DOUBLE  PRECISION  S< 12) » XIC» DETC'XNl »XN2»XN3 

0005  DOUBLE  PRECISION  AO 

0004  DATA  YY/1HY/ 

0007  ND1C  ■  0 

0008  900  TYPE  400 

6009  400  FORMAT < '  INPUT  DATA  FILE  NAME  '.*) 

0010  ACCEPT  401 'FILE 

0011  401  FORMAT <10A4> 

0012  CALL  ASSIGN (3 'FILE tO) 

0013  CALL  FDBSET< 3» 'READONLY') 

0014  READ  <  3 » ERR-900 ) DUM 

0015  REWIND  3 

0014  IF<NDIC.EQ.O.  )  GO  TO  2002 

0018  TYPE  432 

0019  ACCEPT  401 'YB 

0020  IF(YB.NE.YY)  GO  TO  433 

0022  432  FORMAT < '  NEW  CHEST  DIMENSIONS  ?''*> 

0023  2002  TYPE  402 

0024  402  FORMAT ( '  APC 'LATC' APAB 'LATAB  =?  '.•> 

0025  ACCEPT  403 'APC 'XLC 'APAB 'XLAB 

0024  TYPE  419' APC » XLC ' APAB 'XLAB 

0027  419  FORMAT < '  APC  -'»F4.1''  LC  «''F4.1»'  APAB  ='»F4,1»'  LAB  »' 

l'F4«l> 

0028  NDIC  -  1 

0029  403  FORMAT  < 4F8 « 0 ) 

0030  433  TYPE  500 

0031  500  FORMAT < '  FRC  VOL  ■  T  ''*) 

0032  ACCEPT  501  'VFRC 

0033  501  FORMAT <F8. 3) 

0034  TYPE  502  »VFRC 

0035  502  FORMAT ('  FRC  VOLUME  IN  L  ■  '»FB.3> 

0034  IC«1 

0037  DO  54  I  -  1 ' 12 

0038  54  S(I>  -  0 

0039  AO-  APCBXLC 

0040  PI-3. 1415924 

0041  40  READ(3'END«50)  < VVV'D1CC'D2CC'D1AC»D2AC'DUM) 

0042  V  -  VW 

0043  D1C  -  D1CC 

0044  D2C  -  D2CC 

0045  D1A-D1AC 

0044  D2A-D2AC 

0047  D1C-D1C  4  APC 
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0048 

D2C*  D2C  +  XLC 

0049 

D1A  -  D1A  +  APAB 

0050 

D2A  *  D2A  +  XLAB 

0051 

DC  *  PIBD1CBD2C/4 « 

0052 

DA  *  PIBD1ABD2A/4. 

0053 

678 

FORMAT ( 2X  r F8 . 3 ) 

0054 

V  «<  V  +VFRC >41000 . 

0055 

S(l)  -  S<1)  4  DCBDC 

0056 

S(2)  -  S(2>  +  DC/DABDC 

0057 

S ( 3 )  =  S ( 3)  +  DC 

0058 

S<4)  *  S(4)  +  V*DC 

0059 

S(6)  •  S(6)  +  DCBDC/DA/DA 

0060 

S<7)  *  S(7)  +  DC/DA 

0061 

S<8)  =  S(8)  +  DCBV/DA 

0062 

S<11)  -  SOI)  +  1 

0063 

S<12)  ■  S( 12)  +  V 

0064 

IC  *  IC  +  1 

0065 

GO  TO  60 

0066 

50 

S<5)  -  S<2) 

0067 

S<9)  -  S(3) 

0068 

S<10)  -  S<7) 

0069 

TYPE  99f<SO)fI  =  1f12> 

0070 

DETC  -  S0)B(S<6)BS01)  -  S<7)*S<10)) 

0071 

DETC  -  DETC  +  S<9)B<S<2)BS<7)-S(3)BS<6> ) 

0072 

DETC  -  DETC  -S<5)B<S<2)BS01 >-S(3)*S(10> 

0073 

XN1  -  S(4)*(S(6)*SOl)-S(7)*SO0)) 

0074 

XN1  -  XN1  +  S(12)I(S(2)*S(7)-S(3)*S(6)) 

0075 

XN1  -  XN1  -  S(8)*(S(2)*S(11)-S(3)*S(10) ) 

0076 

XN2  «  8(1)*(S(7)*S(12)-S(8)*S01)  ) 

0077 

XN2  ■  XN2  +  S(9)*<S(3)«S(8)-S(4)*S(7)) 

0078 

XN2  ■  XN2  -  S<5)*<S(3)*S(12)-S(4)*S01)) 

0079 

XN3  ■  8(l)t(S(8)*S(10)-S(6)»8(12)) 

0080 

XN3  -  XN3  +  S<9)*<S(4)*S(6)-S(2)»S<8)> 

0081 

XN3  -  XN3  -  S(5)*(S(4)*S<10)-S(2)*SO2>> 

0082 

IC  -  IC  -1 

0083 

TYPE  99'DETCfXN1fXN2'XN3 

0084 

99 

FORMAT (2Xf4E15«8) 

0085 

HO  >  XN1/DETC 

0086 

VA  -  XN2/DETC 

0087 

VC-  XN3/DETC 

0088 

TYPE  71 

0089 

71 

FORMAT <'  V  -  HOBAC  -  VABAC/AA  -  VC  ' > 

0090 

TYPE  70fH0fVAfVC 

0091 

70 

FORMAT <  *  HO  'E15.8f'  VA  'fE15.8f'  VC 

0092 

TYPE  1000 

0093 

1000 

FORMAT < '  PRINT  EVERY  ?TH  PTT'fB) 

0094 

ACCEPT  1001 fNPO 

0095 

1001 

FORMAT < 13) 

0096 

IC  -  1 

0097 

SER  -  0 

0098 

8X  -  0 

0099 

8X2  -  0 

0100 

8Y  -  0 

0101 

8Y2  -  0  _34_ 

' fE15.8) 
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0102 

0103 

0104 

0105 

0106 

0107 

0108 

0109 

0110 

0111 

0112 

0113 

0114 

0115 

0116 

0117 

0118 

0119 

0120 

0121 

0122 

0123 

0124 

0125 

0126 

0127 

0128 

0130 

0131 

0132 

0133 

0134 

0135 

0136 

0137 

0138 

0139 

0140 

0141 

0142 

0143 

0144 

0145 

0146 

0147 

0148 

0149 

0150 

0151 

0152 


SXY  -  0 
REWIND  3 
TYPE  490 

490  FORMAT  <  // '  PO I NT ' »  5X  » ' MEASURED ' »  5X  » ' CALC ' 1 8X  f ' ERR ' ) 

424  READ  (  3  f  END-430  > ( VW f  D1 CC  f  D2CC  f  D 1 AC  f D2AC » DUM ) 

D1C  -  D1CC 

D2C  -  D2CC 

D1A  -  D1AC 

D2A  »  D2AC 

VM  -  VW 

D1C-D1C4APC 

D2C-D2C4XLC 

D2A-D2A4XLAB 

D1A-D1A4APAB 

DC  -  PI*DlC*D2C/4. 

DA-  PI*DlA*D2A/4. 

VV  -  HO* (DC)  -  VA*< DC/DA) -VC 
VM  -  1000. *<VM  +  VFRC) 

ER  -  VV  -  VM 
SX  -  SX  +  VM 
8X2  «  SX2  4  VM*VM 
SY  »  SY  4  VV 
SY2  -  SY2  4  VV*VV 
SXY  -  SXY  4  VM*VV 
NN  -IC  4  NPO-1 

qcd  .  CCD  x  rp*ro 

IF ( NN/NPOBNPO . NE . NN ) GO  TO  1003 
TYPE  425 f  ICfVMfVVfER 

425  F0RMAT(3X» I4*5X» F8. 1 f5XfF8. 1 f5XfF8. 1 > 

1003  IC  -  IC  4  1 

GO  TO  424 

430  TYPE  427 fSERfSER/< IC-1 ) f  <SER/( IC-f)  )**.5 
IC  -  IC  -  1 
SXY  •  SXY  -  SX*SY/IC 
SXY  -  8XY/< IC-1 ) 

SX  -  SX2  -  SXtSX/IC 
SX  -  SX/<IC-1> 

SX  -  SXM.5 

SY  -  SY2  -  SY*SY/IC 

SY  -  SY/< IC-1 ) 

SY  -  8Y**.5 
CC  -  SXY/8X/SY 
TYPE  IIOOfCC 

1100  FORMAT < '  CORRELATION  COEFFICIENT  'fF8.4> 

427  FORMAT < / '  ERROR  80  ' »E12.4»3X» 'ERRSQ/N  'fE12.4»'  RMS 

1E12.4) 

400  FORMAT < '  ANY  MORE  T'»*> 

TYPE  400 

CALL  CLOSE <3> 

401  FORMAT (Al) 

ACCEPT  401 »YB 


0153  XF(YB.EQ. YY)GO  TO  900 

0155  STOP 


FORTRAN  IV  V01C-03C 

C0RE=19K,  UICaC 160>  1 1 


MON  20-AUG-79  09:15528  PAGF  001 

CALC  »CALC/-SP=CALC 


C********************************************************* ******** 

c 

c 

C  PROGRAM  CALC 

C 

C  PREPARES  MODEL  FITS  FOR  VERSATEC  PLOTTING 

C 

C 

C*********************************************************** ******** 

0001  DOUBLE  PRECISION  DC  f  DA  >  VV  t HO » HO  t VA » VC 

0002  DOUBLE  PRECISION  VMtAPC»XLC,APAB»XLAB»VFRC.Pl>CONS»SER>SX 

0003  DOUBLE  PRECISION  SX2f SY»SY2»SXY»DCO»DAO»DlCOrDlAO»D2CO»D2AO 

0004  DOUBLE  PRECISION  D1C*D2C»D1A»D2A»D2D,ER,CC 

0005  DIMENSION  VV<999) »VM<999> 

0006  DIMENSION  FILE< 10) 

0007  DOUBLE  PRECISION  S< 12) »XIC»DETCrXNl »XN2.XN3 

0008  DOUBLE  PRECISION  AO 

0009  DATA  YY/1HY/ 

0010  TYPE  4000 

0011  4000  FORMAT < '  DO  YOU  NEED  MODEL  DESCRIPTION?  ',*> 

0012  ACCEPT  401 tY 

0013  IF(Y.NE.YY)  GO  TO  4005 

0015  TYPE  4010 

0016  4010  FORMAT <////'  MODEL  CALCULATIONS  FOR  ELIPTIC  CYLINDERS' 

1////) 

0017  TYPE  4006 

0018  4006  FORMAT ( '  TYPE  1  V  =  A*<AC-ACO>  -  B*( AAB-AABO)  '/ 

1'  TYPE  2  V  *  A* ( D1C < DI AM ) -D1C0 )  -  B*(D1A(DIAM)  -  D1A0)'/ 

2'  TYPE  3  V  *  A  -  B<<APAB<DI AM)-  APABO)'/ 

3'  TYPE  4  V  -  A* < APC < DI AM ) -APCO )  -  B'/ 

4'  TYPE  5  V  ■  A* ( AC )  -  B*<AC/AAB>  -  C'/ 

5'  TYPE  6  V  ■  A*<AC-ACO>  -  B*<*AC/AAB-AC0/AAB0>  -  C'/ 

6'  TYPE  7  V  -  A*AC  -  B*AAB  -  C'///> 


0019 

4005 

CONTINUE 

0020 

NDIC  -  0 

0021 

900 

TYPE  600 

0022 

600 

FORMAT <'  INPUT  DATA  FILE  NAME  '»*) 

0023 

ACCEPT  601 »FILE 

0024 

601 

FORMAT  < 10A4) 

0025 

CALL  ASSIGN! 3»FILErO> 

0026 

CALL  FDBSET <  3 1 ' READONLY ' ) 

0027 

READ  <  3  f ERR-900 ) DUM 

0028 

4999 

REWIND  3 

0029 

TYPE  6000 

0030 

6000 

FORMAT < '  MODEL  TYPE?  '»*) 

0031 

ACCEPT  1001 »MTY 

0032 

TYPE  6001 »  MTY 

0033 

6001 

FORMAT <//'  COMPARISON  FOR  MODEL  TYPE'»I3//> 

0034 

TYPE  9876 

0035 

9876 

FORMAT < '  PLOT  COMPARISON  ?  '»♦> 

0036 

ACCEPT  401 f YYYY 

0037 

IF< YYYY.NE* YY)GO  TO  9875 

0039 

TYPE  9874 

0040 

9874 

FORMAT < '  FILE  NAME?  '*$) 

-37- 


MON  20-AUG-79  09 J 15 528  PAGE  002 

CALC » CALC/-SP=CALC 


FORTRAN  IV  V01C-03C 

CORE-1 9K»  UIC-C 160» 1 J 


0041 

ACCEPT  601 r FILE 

0042 

CALL  ASSIGN <4*FILE»0) 

0043 

CALL  FDBSET  !  4  #  ' NEW ' ) 

0044 

9875 

CONTINUE 

0045 

IF!NDIC.EQ.O.  )  GO  TO  2002 

0047 

TYPE  432 

0048 

ACCEPT  401 »YB 

0049 

IF! YB.NE. YY)  GO  TO  433 

0051 

432 

FORMAT < '  NEW  CHEST  DIMENSIONS  ?'»♦) 

0052 

2002 

TYPE  602 

0053 

602 

FORMAT !'  APC  »  LATC  *  APAB  »  LAT AB  =?  '»*> 

0054 

ACCEPT  603 • APC  r XLC  r APAB  » XLAB 

0055 

TYPE  6 1 9  >  APC • XLC  »  APAB  r  XLAB 

0056 

619 

FORMAT ! '  APC  ='»F6.1»'  LC  ='»F6.1>'  APAB  ='»F6.1» 
1 »F6. 1 ) 

0057 

NDIC  -  1 

0058 

603 

FORMAT (4F8.0) 

0059 

433 

TYPE  500 

0060 

500 

FORMAT ( '  FRC  VOL  -  T  '»•> 

0061 

ACCEPT  501  * VFRC 

0062 

SOI 

FORMAT <F8. 3) 

0063 

TYPE  502  »VFRC 

0064 

502 

FORMAT <'  FRC  VOLUME  IN  L  »  'rF8.3> 

0065 

IC=*1 

0066 

DO  56  I  =  In  12 

0067 

56 

S<I)  *  0 

0068 

AO-  APC4XLC 

0069 

PI-3.1415926 

0070 

GO  TO  1234 

0071 

678 

FORMAT (2X.F8. 3) 

0072 

1300 

FORMAT (3<2X*E12»5)) 

0073 

1234 

CONTINUE 

0074 

TYPE  1235 

0075 

1235 

FORMAT < '  ArBfC  *?  ' ?*) 

0076 

ACCEPT  1236>H0fVA»VC 

0077 

1236 

FORMAT! 3F12.0) 

0078 

CONS  -  0 

0079 

999 

TYPE  70*H0t VA»VC 

0080 

70 

FORMAT  <  *  A  'E15.8* '  B  '.E15.8.'  C  '»E15.8) 

0081 

76 

FORMAT! F8.0) 

0082 

TYPE  1000 

0083 

1000 

FORMAT!'  PRINT  EVERY  TTH  PTT',9) 

0084 

ACCEPT  1001 »NPO 

0085 

1001 

FORMAT! 13) 

0086 

IC  -  1 

0087 

8ER  -  0 

0088 

SX  -  0 

0089 

8X2  -  0 

0090 

8Y  ■  0 

0091 

8Y2  ■  0 

0092 

8XY  «  0 

0093 

DCO  -  PI*APC*XLC/4. 

0094 

DAO  ■  PX8APAB8XLAB/4. 

0099 

D1C0  •  APC  ”38- 

'  LAB  =  ' 


■38 


FORTRAN  10  V01C-03C 

C0RE-19K*  UIC-C160rl3 


HON  20-AUG-79  09:15:28  PAGE  003 

CALC , CALC/-SP-CALC 


0096  D1A0  -  APAB 

0097  D2C0  -  XLC 

0098  D2A0  =  XLAB 

0099  REWIND  3 

0100  TYPE  490 

0101  490  FORHAT<//'  POINT' »5Xt 'MEASURED' *5X» 'CALC' »8X» 'ERR' > 

0102  424  READ  <  3  *  END-430 ) ( VVM *  DD1 »  DD2  *  DD3  *  DD4  *  DUM ) 

0103  D1C«  DD1 

0104  D2C-DD2 

0105  D1A-DD3 

0106  D2A-DD4 

0107  VM(IC)  »  WM 

0108  D1C-D1C+APC 

0109  D2C-D2C+XLC 

0110  D2A-D2A+XLAB 

0111  D1A-D1A+APAB 

0112  DC  -  PI«DlC*D2C/4. 

0113  DA-  PI*DlA*D2A/4. 

0114  IF<MTY.E0.1>  VV<IC>  «  HO*(DC  -  DCO)  -  VAt ( DA-DAO ) +  CONS 

0116  IF<HTY«EQ*2>  VV(IC)  -  H0*(D1C-D1C0)  -  VAt(DlA-DlAO>+  CONS 

0118  IF(MTY.EQ.3>  VV(IC>  -  HO  -  VA*(D1A-D1A0> 

0120  IF(MTY.E0.4>  VV<IC>  «  HOt(DlC-DlCO)  -  VA 

0122  IF(HTY.EO.S)  VV(IC>  -  HOtDC  -  VAtDC/DA  -  VC 

0124  IF(MTY.E0.6>  VV(IC)  -  HO*(DC-DCO)  -VA*(DC/DA-DC0/DA0>  -  VC 

0126  IF(MTY.E0.7)  VV(IC)  -  HOtDC  -  VAtDA  -  VC 

0128  VM(IC)  «  1000.t(VM(IC)  4  VFRC) 

0129  ER  «  VV(IC)  -  VM(IC) 

0130  8X  -  SX  +  VM(IC) 

0131  8X2  -  8X2  4  VM< IC)tVM< IC) 

0132  8Y  «  SY  +  VV<IC) 

0133  SY2  -  8Y2  4  VV< IC)*VV( IC) 

0134  8XY  -  SXY  +  VM<IC)*VV<IC) 

0135  NN  -1C  4  NPO-1 

0136  8ER  -  8ER  4  ERtER 

0137  IF  <  NN/NPOtNPO . NE . NN ) GO  TO  1003 

0139  TYPE  425*  ICrVM(lC) »VV(IC>  »ER 

0140  425  FORMAT (3XfI4f5X*F8«lf5X?F8«l*5XrF8»l> 

0141  1003  IC  -  IC  4-  1 

0142  00  TO  424 

0143  430  TYPE  427*SER»SER/< IC-1 ) * <SER/( IC-1 ) >tt.5 

0144  IC  -  IC  -  1 

0145  C0N8  -  <8X-8Y)/IC 

0146  TYPE  9000 

0147  9000  FORMAT:'  0FF8ET  T  '»•) 

0148  ACCEPT  401 »HV 

0149  IF(HY.NE.YY)  00  TO  9001 

0151  IF(HTY*EQ«3>  HO  -  HO  +  CONS 

0153  IF<MTY«EQ.4>  VA  -  VA  -C0N8 

0155  IF<MTY'EQ»5.0R'HTY'EQ.6'0R'MTY'EQ.7>  VC  -  VC-C0N8 

0157  TYPE  77*C0N8 

0158  77  FORHATC'  NEW  0FF8ET  ' *F12.3>3XfF12.3) 

0159  9001  SXY  -  SXY  -  8X88Y/IC 

0160  SXY  -  SXY/ (IC-1) 

0161  SX  -  8X2  -  SXSSX/IC 


■39- 


FORTRAN  IV 

V01C-03C  HON  20-AUG-79  09115528  PAGE  004 

CORE- 

19Kt  UIC-C160.1D  CALC » CALC/-SP-CALC 

0162 

SX  -  SX/<IC-1> 

0163 

SX  -  8X4*. 5 

0164 

SY  -  8Y2  -  SY*SY/IC 

0163 

SY  -  SY/< IC-1 ) 

0166 

8Y  -  8Y**«3 

0167 

CC  •'  8XY/8X/SY 

0168 

TYPE  1100.CC 

0169 

TYPE  1200 

0170 

1200 

FORMAT < '  REPRINT  DATA  WITH  NEW  CT  ',*> 

0171 

ACCEPT  401 »YB 

0172 

IF(YB.EO.YY)  00  TO  999 

0174 

1100 

FORMAT < '  CORRELATION  COEFFICIENT  '.FB.4> 

0175 

427 

FORMAT ( / '  ERROR  SQ  ' .E12.4.3X. 'ERRSO/N  '.E12.4.'  RMS 
1E12.4) 

0176 

400 

FORMAT ('  ANY  MORE  ?'  .*> 

0177 

IF(YYYY.NE.YY)  00  TO  8755 

0179 

100-0 

0180 

TYPE  8760 

0181 

8760 

FORMAT < '  EVERY  ?TH  POINT  PLOTTED?  '»*> 

0182 

ACCEPT  1001  »NP 

0183 

NTOT  -  IC/NP+  1 

0184 

IF <  <NT0T-1  )*NP.EQ.  IONTOT  -  NTOT-1 

0186 

WRI TE < 4 . 8765 ) NTOT , 100 

0187 

100  -  4 

0188 

8763 

FORMAT < 13.'. '.Il» ') 

0189 

DO  8764  I  «  l.IC.NP 

0190 

XIC  -  I 

0191 

WRITE<4»8763)XIC. VV< I ) 

0192 

8764 

CONTINUE 

0193 

8763 

FORMAT (2X.F8.0. ' . 'F9.4> 

0194 

WRITE (4. 8765) NTOT .100 

0193 

DO  8762  I  -l.IC.NP 

0196 

XIC  -  I 

0197 

8762 

WRITE (4. 8763) XIC. VM( I ) 

0198 

CALL  CL08E  <  4  > 

0199 

8733 

TYPE  5000 

0200 

3000 

FORMAT < '  ANY  MORE  ON  THIS  FILE?  '.*> 

0201 

ACCEPT  401. Y 

0202 

IF<Y.EO« YY)00  TO  4999 

0204 

CALL  CL08E(3) 

0203 

TYPE  400 

0206 

401 

FORMAT (Al) 

0207 

ACCEPT  401. YB 

0208 

IF(YB.EQ.YY)00  TO  900 

0210 

STOP 

0211 

END 

MWO  M44M 


-40. 


